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Abstract We solve the Ginzburg-Landau equation (GLE) for the mesoscopic supercon- 
ducting thin film of the square shape in the magnetic field for the wide range of 
the Ginzburg-Landau parameter 0.05 < K e ff < oo. We found that the phase 
with the antivortex exists in the broad range of parameters. When the coherence 
length decreases the topological phase transition to the phase with the same total 
vorticity and a reduced symmetry takes place. The giant vortex with the vorticity 
m = 3 is found to be unstable for any field, £/a and K e ff > 0.1. Reduction of 
K e f f does not make the phase with antivortex more stable contrary to the case 
of the cylindric sample of the type I superconductor. 



Introduction 

Advances in nanotechnology and constantly shrinking semiconductor devices 
have motivated researches to study properties of mesoscopic superconducting 
samples. One line of research in this field has focused on the problem of the 
phase transitions in the mesoscopic superconducting sample under the influ- 
ence of the external magnetic field. [ 1 ] We will focus on the case when the 
size of the sample a ~ £, A, with £ and A being the superconducting coherence 
length and the London penetration depth respectively. In that case there are 
only a few vortices in the sample. The standard Abrikosov approach [ 2 1 must 
be modified because of the strong influence of sample boundaries. Thermo- 
dynamics of the system is determined by the short-range repulsion of vortices 
and interaction of vortices with boundaries. 
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Recently it was shown that the influence of boundaries can lead to stabi- 
lization of the vortex-antivortex molecules in mesoscopic samples. [1] Anal- 
ysis of the linearized Ginzburg-Landau equation (GLE) has shown that such 
molecules appear at particular values of the external magnetic field depending 
on the sample shape and size[3|. The solution of the GLE in the limit of the 
extreme type-II superconductor shows that such molecules have a very shallow 
minimum in the free energy[4 ( 5| and are very sensitive to the change of the 
sample shape (6). 

In a square mesoscopic thin film with the total vorticity m = 3 the sym- 
metric solution with four vortices and one antivortex is the solution of the 
linearized GLE with the lowest free energy! 1 1- According to ref. |@), away 
from the H C 2 line the giant vortex with vorticity m = 3 is stable and has the 
lowest free energy. This implies that a topological phase transition without 
change of the vorticity and without a reduction of the symmetry should take 
place with the change of the external field or/and the coherence length away 
from the critical-field line. 

Here we report the results of the extensive studies of different kinds of phase 
transitions for the thin film of the square shape. We focused to the region 
4 < a/£ < 8 and K e ff > 0.05 (see also[7|). Here K e jf = A 2 /(f£ where d is 
the film thickness. We found that the antivortex phase with m = 3 is stable 
in a broad range of parameters. The region of stability of the phase does not 
depend strongly on the value of the parameter K e //. The energy gain due to 
the antivortex formation is much smaller then the energy difference between 
two phases with different vorticities. The giant vortex with m = 3 is unstable 
for any field, £/a and K e fj > 0.1. Phase transition to the phase with three 
separated vortices takes place when £/a is driven away from the critical field 
line. The reduction of n e f f does not stabilize the antivortex phase for the thin 
film sample in the contrast to the case of the cylindric sample ref. [8 ]. 

1. Formalism and Solution 

GLE for the normalized complex order parameter tp has the following form: 


9-77- A 

£ 2 (*v + ^— ) 2 ^-V + V#l 2 = o (i) 

here <E>o is the flux quantum, A is the vector potential and H = V x A the mag- 
netic field. We split the vector potential into the external part due to external 
currents, A ext , and the internal part due to the response of the superconducting 
film, Ai nt . The second GLE equation for the total vector potential reads: 

V x V x A = - A(f - W ) - (2) 
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In addition to Eq.(l) we assume the boundary condition for the superconductor 
insulator junction on the sample edges: 

(iV + ^)-n^ = 0, (3) 

where n is the vector normal to the surface of the sample. 

As it was described in ref.f?) we introduce N x N discrete points on the 
square and rewrite Eq. (1) in the form of the nonlinear discrete Schredinger 
equation: 

^ i+M V>i+i - e(i)i M V>i - + V'ilV'il 2 = 0, (4) 
l 

where the summation index 1 = (±1,0), (0, ±1) points toward the nearest 
neighbors and t hti = (£N/a) 2 exp(i(f> hti ) and iX)i = -|jJJiA(r)(ir. 
The boundary conditions are included in the discrete nonlinear Schredinger 
equation as in ref.[4| where tpi = if i is outside of the sample and e(i) = 
4 - S ix> i - S ixjN - 5 iyt i - 5 iytN where i = (i x = 1, . . . , N, i y = 1, . . . , N). 

After discretization of Eq. (2) we obtain the exact expression for the internal 
part of the vector potential 171: 

n 

where 

Ji = 47rA ° Q jV ^( ex PH0i44v,i)#V'i+lv - 

exp(-i0 i _i v>i )^iV'i-i v ) (6) 
with v G {x, y} and l x = (1, 0), l y = (0, 1) and 

K(n) = r dxdy - ™{n x x) C o S {n y y) 

V ; 2vr 2 ai - 2cos(x) - 2cos(y) 

Here we should point out that Eq.(5) contains 2D integration only[7|. All 
dependence on the thickness of the sample appears through the parameter 
X e ff = X 2 /d (d is the thickness of the sample)! 9 1. This is important dif- 
ference from the case of the cylindric sample considered in the Ref.[8 1 where 
the function K (n) is essentially different. 

The numerical self consistent solution of the problem is obtained by iterat- 
ing the solution of the nonlinear equation for the order parameter Eq.(4) and 
calculations of the current and the vector potential Eqs.(5,6). We used two 
ways of solving Eq.(4). The first is similar to that reported in ref.[4| and cor- 
responds to the iterative solution of the linearized Eq.(4). The second relies on 
the fact that Eq.(4) represents the Euler equation for the free-energy functional 
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Figure 1. The modulus of the order parameter at different magnetic fields. The parameter 
<E> is the total external magnetic flux through the sample. The inset in d) shows the central region 
in an expanded scale. Note the logarithmic intensity scale. 



with included boundary conditions. Eq.(4) was therefore solved by the di- 
rect minimization of the corresponding functional using the conjugate- gradient 
method. Both techniques gave identical results. 

2. Results 

We investigate first the phase diagram in the regions where the solution with 
one antivortex and four vortices and one antivortex and eight vortices has been 
reported. In Fig. 1 we present the change of a spatial pattern of the modulus 
of the order parameter when the vorticity changes from m = 2 (Fig. la) to 
m = 4 — 1 (Fig. lb) and from m = 6 (Fig.lc) to m = 8 — 1 (Fig. Id). 

In the cases m = 2 and 6 (Fig. la and Fig. lb) we observe a giant double 
vortex in the centre. In the cases with m = 41 — 1 with I = 1, 2 instead of 
the giant triple vortex the symmetry induced square pattern of four vortices 
with the antivortex in the middle forms. The results of the calculations show 
that the region of the phase diagram where the symmetry induced antivortex 
solution has the lowest energy is broader than expected from the solution of 
the linearized GLE. As it is shown in Fig. 2 for K e f f = oo the antivortex 
phase is stable up to (a/£) 2 ~ 55, depending on ^/^o- F° r a finite n e ff this 
region shifts to the higher field as (a/£) 2 increases (see Fig.2). In the Fig.3 
we presented the calculated value of the sample free energy as a function of 
external field for (a/£) 2 = 35 and for two values of n e ff = 1, oo. Topological 
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Figure 2. The calculated phase diagram. Different phases are marked with icons schemati- 
cally indicating the vortex pattern where the full dot represents a vortex, the open dot reperesents 
an antivortex and the larger full dot represents a double vortex. The full symbols and continous 
lines represent the phase boundaries for n e ff = oo while the open symbols and dotted lines 
reperesent the phase boundaries for n e ff — 1. In the later case only the phase boundaries of 
the region with the total vorticity 3 are shown. 
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Figure 3. The free energy as a function of magnetic field for different vorticities m. 
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phase transitions with the changes of the total vorticity Am = 1 are clearly 
seen. Reduction of n e j t leads to the shift of the transition point to the higher 
field. 




Figure 4. The magnitude of the order parameter around the transition where one vortex 
anihilates with the antivortex as a function of (a/£) 2 at the constant magnetic field. 

The interesting behavior is observed when the external field is fixed and 
(a/£) 2 increases. Close to the H C 2 the lowest minimum of the free energy 
corresponds to the solution with the vorticity m = 4 — 1 with the antivortex 
in the center of the square. Present calculations do not confirm the existence 
of the giant-vortex solution with m = 3 in this region of the phase diagram 
as reported previously |4]. The difference is due to increase of the number 
of discrete points N enabling detection of the antivortex. With increase of 
(a/£) 2 away from the H C 2 line the phase transition to the multivortex state 
with the same vorticity (m = 3) and a lower symmetry takes place (see Fig. 
2). In general, the free energy depends on the vorticity m = n + — ra_ and 
the total number of vortices in the system n = n + + n_. The transition at 
(°/£) 2 ~ 55 and ^/^o ~ 5.5 takes place at the constant vorticity m = 3 with 
the change of n from 5 to 3. The transition is therefore not only characterized 
by an order parameter, but also by the change of the number of vortices at 
the constant total vorticity m, suggesting that the transition is close to the first 
order. This statement is confirmed by the observation that above the transition 
point, (a/£) > {a/£) cr it, both solutions with m = 3 and m = 4 — 1 coexists. 
Since near the transition the free-energy difference between the phases with 
the same vorticity m and different n is small it is difficult to determine the 
phase boundary between phases with m = 4 — 1 and m = 3 accurately. The 



Phase Transitions in Mesoscopic Superconducting Films 



7 



transition could be easier observed by calculating the two component order 
parameter r/ x = J x\ifj(x,y)\ 2 dxdy, r\ y = J y\^(x,y)\ 2 dxdy shown in Fig.4. 




Figure 5. The vortex- antivortex distance a) and the magnetic moment of the sample b) as 
functions of the parameter l//t e //. In (a) error bars represent the grid spacing and the solid line 
the exponential fit discussed in the text. 

Close to the H C 2 line the repulsion of vortices from the boundaries and at- 
traction of the 4 vortices to the antivortex stabilizes the phase with m = 4 — 1 
and small vortex-vortex distances. At smaller value of £ the repulsion from 
the boundaries decreases and one vortex annihilates with the antivortex. As 
a result, the repulsion between the remaining vortices increases leading to an 
increase of the order parameter with a further decrease of £. 

Increasing the external field up to 3>/<]?o = H-6 leads to the stabilization of 
the phase with total vorticity m = 7. Near the H C 2 line similar to the case with 
m = 4 — 1 the solution with m = 8 — 1 is realized (Fig. Id). When (a/£) 2 
increases in a complete analogy to the case with m = 4 — 1 the phase transition 
to the phase with m = 7 and a similar order parameter takes place. Here also, 
both solutions with m = 8 — 1 and m = 7 coexists above (a/^) cr u indicating 
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that the transition is close to the first order. We believe that the situation is 
quite general for the case of arbitrary m = 41 — 1 for I = 1, 2, 3, .... 

At the end we would like to discuss the dependence of the stability of the 
antivortex phase at small ft e // . According to the arguments of Ref . [ 8 1 , at small 
k the vortex-vortex interaction changes the sign making the antivortex phase 
more stable. As a result, the average distance between vortices in the middle 
of the square increases as well. In order to verify this conjecture for the thin 
film sample we plot in Fig. 5 the vortex-antivortex distance tq as a function of 
l/n e ff. The distance decreases with the decreasing K e ff- For n e fj < 0.1 the 
distance is smaller than the grid spacing a/N so we can not resolve separate 
vortices. We find that tq oc exp (— A/\ e ff) with A ~ a. The situation is just 
opposite to that reported in ref.[ 8 1. We believe that in the case of the thin film 
of the square shape the reduction of k does not stabilize the phase with the 
antivortex. 




Figure 6. The magnetic field in the film in the case of (a) giant vortex with m — 2, (b) the 
antivortex solution with m = 4 — 1 and (c) three separate vortices with m — 3. Here Hq is the 
external magnetic field. 
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It is interesting to note differences between samples of different shapes. For 
the cylindric shape the giant vortex phase with any vorticity is always stable 
close to the H C 2 line| 10 11]. According to Ref.[ 8 1 for the mesoscopic triangle 
the giant vortex state with m = 2 is metastable and the solution with the 
antivortex (m = 3 — 1) is stable. For the case of the square shape the giant 
vortex solution with m = 3 is never stable for K e ff > 0.1. For K e ff < 0.1 the 
limited grid spatial resolution prevented us to distinguish the solution with the 
antivortex from the possibly (meta)stable giant-vortex solution. 

Finally, let us discuss the possibility to detect the state with the antivortex 
experimentally. Calculation of the magnetic field in the sample shows that the 
magnetic field has a local minimum in the center of the sample also for the 
giant vortex solution with m = 2. The local minimum observed for the an- 
tivortex state with m = 4 — 1 is therefore not due to the antivortex formation 
(Fig. 6) but due to a particular distribution of the current in the sample. There- 
fore, imaging of the magnetic field distribution cannot provide an evidence for 
the antivortex. The magnetic field for the multivortex solution with m = 3 has 
3 well separated maxima that break the four fold rotational symmetry of the 
sample allowing a direct imaging of vortices. Since the antivortex state cannot 
be detected directly the observation of a hysteresis in the vicinity of the transi- 
tion line from the m = 4 — 1 antivortex state to the m = 3 multivortex state 
could suggest that the symmetric phase is indeed the phase with the antivortex. 
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